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Cross-correlation of 2MASS and WMAP3: Implications for 
the Integrated Sachs- Wolfe effect 
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ABSTRACT 

We perform a cross-correlation of the Cosmic Microwave Background (CMB) using 
the third year Wilkinson Microwave Anisotropy Probe (WMAP) data with the 2 
Micron All Sky Survey (2MASS) galaxy map (about 828 000 galaxies with median 
redshift z 0.07). One motivation is to detect the Integrated Sachs- Wolfe (ISW) effect, 
expected if the cosmic gravitational potential is time dependent; for example, as it is 
in a flat universe with a Dark Energy component. The measured spherical harmonic 
cross-correlation signal favours the ISW signal expected in the concordance ACDM 
model over that of zero correlation, although both are consistent with the data within 
2cr. Within a flat ACDM model we find a best fit value of fi A = 0.85 and il A < 0.89 
(95% CL). The above limits assume a galaxy bias b g (^p^) w 1.40 ± 0.03, which we 
derived directly from the 2MASS auto-correlation. Another goal is to test if previously 
reported anomalies in the WMAP data are related to the galaxy distribution (the so- 
called "Axis of Evil" - AoE). No such anomaly is observed in the 2MASS data nor are 
there any observed AoE correlations between the 2MASS and WMAP3 data. 
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1 INTRODUCTION 

Our currently favoured ACDM cosmological model has 
received on-going confirmation and bolstering over re- 
cent years, especially from recent observations of the cos- 
mic microwave background (CMB) by the Wilkinson Mi- 
crowave Anisotrop y Probe (WMAP) (jSpergel et all 120031 ; 
ISpergel et alj [20061 ). This model postulates that more than 
two thirds of the Universe is composed of 'Dark Energy', 
a mysterious energy with negative pressure. This Dark En- 
ergy has never been directly observed - only inferred. The 
case is compelling because of the very different types of 
observations that require its presence: acceleration of the 
Universe seen by supernovae; joint analysis of the CMB 
with Large Scale Structure (LSS) requiring zero curva- 
ture (|Hinshaw et al.ll2006l ; ISpergel et ail I2006T ) but a mass 
density today Q m ~ 0.25. However, it would still be comfort- 
ing to have a more direct and independent detection of the 
Dark Energy and its effects, especially in light of promising 
alternative theories based on inhomogeneous models which 
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do not require Dark En ergy l|Alnes et alj|2005l : lMoffatJl2006l : 
IVanderveld et al.ll2006l h 

The Integra ted Sachs-Wolfe (ISW) ef- 
fect (|Sachs fc Wolfel Il967r ) provides us with one of the 
cleanest probes of Dark Energy, though it does require 
some assumptions about other cosmological parameters. 
As CMB photons travel through space they pass through 
the gravitational potential wells of LSS. As they fall into 
a well, the photons are blueshifted, and as they climb 
out they are redshifted. In an Einstein-de Sitter universe, 
these effects cancel and no net shift is observed on large 
scales. However, if Dark Energy dominates, large scale 
gravitational potentials decay and there is an overall net 
effect. This secondary anisotropy is called the late-time 
ISW effect. 

On large scales, the ISW effect will add power to the 
CMB anisotropies, by: 



^) =-2 
T J isw 



$'(0?o - v)&,v)dv 



(1) 



where T is the temperature; rj is the conformal time, de- 
fined by dr; = and 770 and 77L are the conformal 
times today and at the surface of last scattering respec- 
tively; n is the unit vector along the line of sight; $(x, rf) 
is the gravitational potential at position x and at confor- 
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mal time r), and $' = Its relative amplitude makes it 
diffic ult to distinguish from the primary anisotropics. How- 
ever, ICrittenden fc Turokl (|l996h proposed using the cross- 
correlation between the LSS and the CMB to detect the 
ISW effect, independently from the intrinsic CMB fluctu- 
ations. In the case where the gravitational potentials de- 
cay, a positive correlation is expected. This means that on 
large scales hot spots in the CMB will correspond to over- 
dense regions in the galaxy distribution. Th is positive cor- 
relation is also expected in open universe s l|Kamionkowskil 
19961 : iKinkhabwala fc Kamionkowskilll999h , whereas a nega- 
tive correlation will occur in closed universes. The Sunyaev- 
Zeldovich effect also produces a negative correlation, but 
this is expected on s maller scales (£ > 20 for z ~ 0.07, see 
lAfshordi et ail (|2004h ). 

In our currently favoured cosmological models, we be- 
lieve the Universe has recently (z < 1) become dominated 
by Dark Energy, which makes the ISW effect a fitting 
probe of our cosmological model, thoug h alternative mod- 
els of gravity pre dict similar signatures (|Carroll et al.ll20o"5l ; 
ISong et al.|[2006h . 

The first cross-correlations between the COBE (Cos- 
mic Microwave Background Explorer) CMB map and 
tracers of the LSS (hard X-ray background and Ra- 
dio s ources) did not reveal a ny s ignificant detec- 
tions l|Boughn fc Crittenden! |2002| . |2003| ); however, there 
have since been a number of reported detections of late- 
time ISW from cross-correlating the first-year WMAP 
data (WMAP1) with: Radio sou rces (|Nolta et al.l |2004 
iBoughn fc Crittenden! |2004 120051^: the h ard X-ray back- 
ground (|Boughn fc Crittenden|[200l.l2005l): the Sloan Digi- 
tal S k y Survey (SPSS) (IScranton et al.ll2003l : iFosalba et all 



I2OO3I ; (Padmanabhan et al.l 120051) ; the 2 Micron All Sky 
Survey (2MASS) (lAfshordi et all |2004| ); the APM Galaxy 
Survey IIFosalba fc G aztafiaga 2004); a combination of the 
above (|Gaztanaga et al.l |2006| V Recently, th e third-year 
WMA P data (WMAP3) was correlated by ICabre et ail 
(2006) with the fourth SDSS data release (DR4), and 
they detected a significan t positive cross-correlation, while 
ICiannantonio et alj |2006l ) cross-correlated it with high red- 
shift SDSS quasars and found a 2a detection. The above 
cross-c orrelations were all per formed in angular or harmonic 
space. iMcEwen et alj (|2006h used a directional spherical 
wavelet analysis and found a positive detection at the 3.9<r 
level. 

A cross-correlation between LSS and the CMB can also 
probe the issue of foreground contamination in the CMB and 
LSS maps. Recent studies have highlighted anomalous fea- 
tures in CMB data that include alignment s of the low-l mul- 
tipoles , and a North-South asymmetry (lLand fc Magueiiol 
l2005al lbl; ICopi et all 120061 ; lEriksen et all |2004| ). Curiously 
these features align with the ecliptic plane, and possibly the 
Supergalactic Plane. Such alignments may, of course, be co- 
incidental, or they may indicate the culprit behind these 
anomalies, such as contamination or some other mor e lo- 
cal secondary effect of t he structure on CMB photons l|Valel 
l2005l ; lRakic et al.ll2006l ). 

2MASS as an infrared survey has the advantage of prob- 
ing the older stellar populations in galaxies and is therefore 
a robust tracer of their mass. Since infrared light penetrates 
more easily through the Galaxy than optical light, 2MASS is 
nearly an all-sky extra-galactic survey. These two qualities 



of 2MASS make it one of the best available tracers of any 
large scale correlations, such as those expected by the ISW 
effect. However, even in the infrared, Galactic extinction and 
stellar contamination will produce a foreground contamina- 
tion in the galaxy distribution which can be correlated with 
emission foregrounds that contaminate WMAP. 

In this paper we perform a harmonic space cross- 
correlation of th e CMB with LSS us i ng the WMAP3 tem - 
perature maps dSpergel et al . 2006; Hinshaw et al. 2006) 
and t he 2MASS galaxy survey ( Jarrett et al . 2000; Jarrett 
|20Q4|). We account for possible correlations from Galactic 
contaminatio n using random si mulat ions. This work updates 
the work of lAfshordi et al. I <|2004h . who cross-correlated 
WMAP1 with 2MASS. We also investigate correlations be- 
tween another statistic, t he so-called "Axis of Evil" (AoE) 
(lLand fc Magueiiol l2005al 1. 

In Section [2] we introduce the 2MASS and WMAP3 
data that we use, and in Section we outline our theory. 
In Section U we discuss our choice of fiducial model. In 
Section [5] we estimate the galaxy bias, 6 g , from the galaxy 
auto-correlation function. In Section[5]we discuss our cross- 
correlation method. In Section[7]we present our results of the 
cross-correlation and discuss their statistical significance. In 
Section [5] we investigate AoE type correlations and in Sec- 
tion [5] we present a discussion of our results. 



2 THE DATA 

In this section we summarize the two data sets used for the 
cross-correlation. 



2.1 The Large-Scale Structure: 2MASS 

We use the publicly available full-sky e xtended source cata - 
log (XSC) of the near-infr a red 2M ASS l| Jarrett et al.ll2000l l. 
Following lAfshordi et all l|2004l ). we divide the galaxies 
into different magnitude bands depending on their K„- 
band isophotal magnitude K20 ( "k_mJ20c" - 20mag/ arcsec 2 
isophotal circular ap. magnitude). We correct these magni- 
tu des for Galact i c exti nction using the IR reddening maps 
of ISchlegel etafl <|l99Sl fT: 

K 20 K 20 - A K (2) 

where A K = 0.367 x E(B - V) ijAfshordi et al.ll2004h . We 
use the same extinction maps to create a mas k that excludes 
region s of sky where the XSC is unreliable. lAfshordi et alj 
(2004) find a limit of A k < 0.05 for which 2MASS is seen 
to be 98% complete for K20 < 13.85. We adopt this method 
herein, and mask areas with Ak > 0.05, leaving 69% of the 
sky with approximately 828,000 galaxies for the analysis. 

For our four K20 magnitude shells, we adopt the 
redshift distribution computed by lAfshordi et al.l (120041 ) 
(parametrized by their Equations 33 and 35). They fit the 
redsh ift distribution from t he 2MASS K20 luminosity func- 
tion (jKochanek et al.ll200ll ) for a three parameter general- 
ized gamma distribution which we use herein. We recall the 
number counts for the four K20 shells and the redshift zo 



1 Corresponding reddening maps can be found 
http: / / astro, berkclcy. cdu /davis /dust /index, html 
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ESTIMATED REDSHIET 
EIT EROM AESHORDI ET AL. 



Figure 1. The redshift distribution for different magnitude cuts 
of 2MASS, estimated from the parametrization in Afshordi et 
al (2004). The red (dotted) line corresponds to the closest shell 
12 < K 2 o < 12.5, the blue (dashed) to the 12.5 < K20 < 13 shell, 
the magenta (dot-dashed) to the 13 < K20 < 13.5 shell, the cyan 
(dot-dot-dashed) to the furthest 13.5 < K20 < 14 shell. The solid 
(black) shows all four shells combined, i.e. 12 < K20 < 14. The 
redshift, zq , at which each distribution peaks and the number 
counts for the four K20 shells are reported in Table [T] 



Table 1 . The number of galaxies for our four K20 shells and for 
all four shells combined. We also recall the peak redshift values 
for the distributions. These values are taken from Afshordi et al. 
(2004). The parametrizations of these distributions are shown in 
Figure [T] 



Magnitude 


Ntot 


zo 


12.0 < K20 < 12.5 


49902 


0.043 


12.5 < K 20 < 13.0 


102947 


0.054 


13.0 < K 20 < 13.5 


217831 


0.067 


13.5 < K 20 < 14.0 


457267 


0.084 


12.0 < K20 < 14.0 


827947 


0.073 



at wh ich the distributions peak, taken from lAfshordi et al.l 
(2004), in Table [1] In Figure [1] we plot the redshift dis- 
tribution for each shell as well as for 12 < K20 < 14. The 
2MASS overdensity field for galaxies with galactic extinction 
Ak < 0.05 is plotted in Figures [2] where it is convolved with 
a gaussian beam of Full Width Half Maximum (FWHM) 
100'. 

We make maps of the 2MASS ove rdensity using the 
Healpix format (|Gorski et al.lll998l . [20051 ), and we use their 
'map2alm' routine to obtain the spherical harmonic coeffi- 
cients. For an incomplete sky, the magnitude of the angular 
power spectrum scales with the survey area and we com- 
pensate for the loss of sky cover by including the factor f s k y 
in Equation 1141 The sky cut will also induce correlations 
between adjacent multipoles. 



2.2 The Cosmic Microwave Background: WMAP3 

We use the third-year data from NASA 's WMAP satellite 
|Hinshaw et al.ll200rj ; [Spergel et al.ll2006l 'Fl. We use the fore- 
ground reduced maps of the Q (41 Ghz),V (61 GHz), and 
W (94 GHz) bands. The foreground reduced sky maps were 
produced by removing a foreground model from the "unre- 
duced" maps. Synchrotron, free-free, and dust emission tem- 
plates were modelled and then subtracted from the single 
year "unreduced" maps. Full three year maps were produced 
by performing a weighted, pixel-by-pixel, mean of the three 
single year maps. This same weighted mean method was 
then used to combine three year maps of the same frequency 
band into a single map for each frequency. Two maps were 
combined to produce the Q and V band maps; four maps 
were combined to produce the W band map. We use the 
Kp2 mask to exclude the Galactic plane and other known 
foreground sources. We also use the WMAP3 ILC map with 
a Kp2 mask. The WMAP3 ILC map convolved with a gaus- 
sian beam with FWHM 100' is plotted in the top part of 
Figure[5]hi Mollweide projection. As with the 2MASS data, 
we use the Healpix routine 'map2alm' to find the spherical 
harmonic coefficients. 

We also return to the WMAP1 data and in Section[U]we 
compare the cross-correlation results for the V band (which 
is the least contaminated band with the best resolution). 
We use the first-year foreground cleaned maps, with the 
Kp2 mask, and inverse-noise coadd the two available V band 
maps. 



3 THEORY 

In this section, we present the formalism for the calculation 
of Auto- and Cross-Correlation Functions (ACF & CCF). 

In what follows, we use r as the comoving distance and 
implicitly as a label of redshift epoch z. They are related for 
a given cosmology by dr = -j^-^dz, where c is the speed of 
light in vacuum. Thus, the growth factor D(z), the growth 
function f(z) (defined below) and the Hubble parameter 
H(z) all have implicit dependences on r. 

For the angular CCF: 



c gT (i) = 4^fo g 



dk 



A 2 (fc) 
k 



Wg(k)W T (k) 



For the angular overdensity ACF: 
A 2 (fc) 



C gg (£) = 4tt^ 
where: 

W g (k) 
W T (k) 
6(r) 



dfc ■ 



k 



W g (k) 



dr @(r)ji(kr)D 



3f2 m ,o-ffo ' ' 
k 2 c 3 

r 2 n c {r) 



dr j e (kr)HD( f - l) 



J dr r 2 n c {r) 
2 Available at http://lambda.gsfc.nasa.gov/ 



(3) 
(4) 

(5) 
(6) 
(7) 
(8) 
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WMAP 




2MASS 




Figure 2. Best viewed in colour. Top: The 3rd year WMAP (WMAP3) Internal Linear Combination (ILC) map, which maps 
temperature anisotropics in the Cosmic Microwave Background. Bottom: The 2 Micron All-Sky Survey (2MASS) overdensity field of 
827,947 galaxies with magnitudes 12 < if 20 < 14. Red shading represents CMB hot spots and 2MASS overdensities and blue shading 
corresponds to CMB cold spots and 2MASS underdensitics respectively. Grey shading corresponds to the K p 2 mask for the CMB and 
for the 2MASS map to galaxies with galactic extinction A k > 0.05. Both maps are produced with Healpix, are in Mollweide projection of 
Galactic Coordinates and are convolved with a gaussian beam of FWHM 100'. The center of the projection represents galactic coordinates 
(I, b) = (0, 0) and the longitudinal coordinate increases leftwards and I = 180° on both left and right edges of the projection. The data 
are publicly available from http://lambda.gsfc.nasa.gov/ (WMAP3, see Hinshaw et al. 2006) and http://www.ipac.caltech.edu/2mass/ 
(2MASS, see Jarrett 2004). 
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Integrals are performed over the wavenumber k (ex- 
pressed in /iMpc -1 ), and the comoving distance r (in 
ft _1 Mpc). The matter power spectrum, P(k), is related to 
the galaxy power spectrum through the linear bias b g which 
we take to be constant on the depth of o ur survey. This is jus- 
tified by the Galaxy conserving model of Ma gliocchetti et alj 
(2000) (their Equation 2) where the bias evolution between 
redshift zero and the mean redshift of 2MASS was found 
to be less than 1%. The selection function, n c (r), for each 
magnitude shell, is shown in Figure [1] 

The linear growth factor, D(z), is given by: 



D(z) oc H(a) 



a(z) 



do' 



{a'H{a')) 3 



(9) 



and is normalized such that D(0) — 1. In Equation [7] the 
redshift at the surface of last scattering is zl ~ 1089. The 
growth function / is given by f(z) = ^J^ xl] a n d can 
be approximated by f(z) ~ 0.t 6 (z) l|Peeblesl (|l993l ). their 
Equation 5.120). The dependence on the Dark Energy den- 
sity and equation of state is negligible at the present epoch 
jLahav et al.lll99ll : IWang fc Steinhardtlll998l ). The Hubble 
constant is parametrized by Ho = 100ft km s" 1 Mpc -1 and 
je(kr) is the spherical bessel function of the 1st kind of order 
I. 

These equations are exact in linear theory, and for t > 
10 we replace the m with the small angl e approximation or 
Limber equation (lAfshordi et al.|[2004 ). which arise from 
the bessel function approximation: 



lim j e (x) 



5 £■ 



21+1 

Equations [3] and [4] then simplify to: 



C gT {£) = 



C gg {£) 



c 3 (£ + 1/2) 2 



a 2 

dr -rr D 2 P 



dr @D 2 H[f-l]P 
'£ + 1/2 



£+1/2 



(10) 

(11) 
(12) 



Reducing the number of integrals significantly reduces com- 
putation time. In Figure [3] we compare the CCF using the 
equation exact in linear theory (Equation [3} with its small 
angle approximation (Equation lll[) . The difference is less 
than 1% from I = 5 upwards. However, for lower multipoles 
the difference is more important; for example, at £ = 2 it 
is of order 10%. The small angle approximation makes the 
assumption that kr ~ (£ + 1/2). From this we see that for 
deeper surveys, the small angle approximation will begin 
to hold at a higher multipole. One should therefore check 
at which multipole the Limber equation begins to hold. 
Throughout this paper we use Equationsl3landl4"l (|lll and ll2l 
when £ > 10) to model our theory. 

In Figure [4] we plot the CCF for different values of 
the matter density, Q m , for flat cosmologies with Dark 
Energy. The correlation increases with increasing We 
note that an ISW effect is also observed in open cosmolo- 
gies without Dark Energy ([K amionkowski & Spcrecll 1 19941 : 
iKamionkowskil 1996; Kinkha bwala fc Kamionkow ski 1999); 
and for low redshifts (z < 2) the signal expected in an 
open cosmology can be grea ter than the s i gnal expected 
in a flat ACDM cosmology (Ka mionkowskil Il996l ). Closed 
cosmologies would produce a negative correlation. In fact, 



Small Angle Approximation 
Exact in Linear Theory 




Figure 3. The expected cross-correlation from the equation exact 
in linear theory (Equation \3\ - solid red line) and using a large- 
angle approximation for the bessel function (Equation llll - dashed 
blue line) . These values are obtained using our fiducial model (flat 
universe with f2 m = 0.30, C„ = 0.05, h = 0.7, as = 0.75) and 
the selection function for galaxies with 12 < K20 < 14. The 
black dot-dashed line is the absolute relative difference between 
the two. The difference is less than 1% from £ = 5 onwards, 
but is considerably larger for lower multipoles. For 1 = 2, for 
example, the difference is of order 10%. The multipole at which 
the small angle approximation begins to hold increases with the 
mean redshift of the survey. 



n=o.2, n.=oi 



Q =0.3, n A =0.7 



n =0.4, n A =o.6 



\ 

s 



Figure 4. The expected ISW signal (fif, = 0.05, h = 0.7, 
erg = 0.75, b = 1.40) for the fourth magnitude bin for flat cosmolo- 
gies with a cosmological constant. Increasing correlation corre- 
sponds to increasing values of f2y\. The dotted line corresponds to 
a cosmology with f2 m = 0.40, the dash-dotted line to f2 m = 0.30 
and the dashed line to O m = 0.20. The ISW effect is strongest 
at scales i < 30. An ISW effect is also observed in open cos- 
mologies without Dark Energy (Kamionkowski & Spergel 1994; 
Kamionkowski 1996: Kinkhabwala & Kamionkowski 1999). 



within the family of cosmological models with arbitrary mat- 
ter density, cosmological constant and curvature, only the 
Einstein-de Sitter universe (tt m = 1, £Ia = 0) will give zero 
correlation. This can be seen from the last term in Equa- 
tion [7] for no correlation, the following condition must be 
satisfied: f(z) ~ Q^(z) = 1 for all redshift. 
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4 CHOICE OF FIDUCIAL MODEL AND 
PRIORS 

In this paper we focus on how to constrain Dark Energy 
from the ISW effect alone. This requires of course assump- 
tions about the other cosmological parameters. To have the 
ISW result independent of the CMB, we prefer to assume 
"round" cosmological parameters which are in accord with 
other cosmological measurements, rather than adopting ex- 
act values from another analysis, e.g., the WMAP3 TT cor- 
relation function. 

Based on inflation, we assume the Universe is flat, there- 
fore Qa = 1 — fim- We also assume the spectral index to 
have the Harrison-Zeldovich value n = 1 (although some 
inflationary models and the recent WMAP3 data suggest 
n w 0.95). Based on the HST key project, we take for 
the Hubble parameter h = 0.7 (the la error bar is about 
10%). For the matter density we take Q m = 0.30, in ac- 
cord with supernovae la data combined with the flatness of 
the universe (we note however th at the recent 2dFGR S data 
(Cole et al. 2005) and WMAP3 i|Spergel et alj|2006ft favour 
fim w 0.25). For the baryon density we assume Qb = 0.05, 
based on Big Bang nucleosynthesis l|Copi et al.lll995l) . 

The normalization of the power spectrum, parametrized 
as as, is still highly uncertain, with reported values in the 
range as ~ 0.75 (e.g . the recent WMAP3 result) to as ~ 1.0 
l|Massev et al.ll2005h . We fix a s = 0.75, and we solve for the 
galaxy biasing 6 g . We note that in linear theory we actually 
constrain the product b s a%, so we can easily scale the result 
for any preferred value of as- 

In principle one should marginalise over the prior as- 
sociated with each parameter. However, as we show in Sec- 
tion [71 the cross-correlation signal is very weak and is actu- 
ally compatible with the null hypothesis of no correlation. 
By widening the priors on 6 g and as, we find this makes 
the signal even less significant. We are aware that widening 
priors for each parameter will reduce the significance of our 
result, but we chose to opt for a more optimistic view, and 
acknowledge that the significance of our result is tentative. 
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Figure 5. The raw (left) and binned (right) 2MASS angular 
power spectrum (solid line) for galaxies 12 < K20 < 14, with 
theory (dashed) for the best fit bias b g = 1.40, and lc gaussian 
error bars (dotted). The expected poisson noise has been sub- 
tracted from the data. We use Q m = 0.30, fit, = 0.05, as = 0.75, 
h = 0.7 as our fiducial model. 



scale-dependence to change much. This issue does not arise 
when considering the CCF, as the ISW effect should only 
arise on linear scales. 

We therefore decide to fit the ACF using a non-linear 
power spectrum for I ^ 50 (to avoid highly non-linear 
scales), and we bin the data into 6 logarithmically spaced 
bins. This reduces correlation between different multipoles 
so that we can assume the bins are independent and have 
gaussian error bars. Hence, scatter about the expected sig- 
nal is just due to cosmic variance, and the likelihood can be 
written: 



- 2 In C = det(M) + d T M _1 d + const 



(13) 



Where the diagonal terms of M contain the variance of C„ 



Whe re the d iagonal tei 
fe.g.: lDodelsonf l|2003h ) 



5 GALAXY BIASING FROM THE GALAXY 
ACF 

We compute the ACF of 2MASS, C gg (£), and use it to con- 
strain the galaxy bias, 6 g , in Equations l4l (and I12|l . Since we 
assume a constant bias across the magnitude shells we make 
the fit using all the galaxies at once, 12 < K20 < 14. The 
redshift distribution for all of 2MASS is plotted in Figure [JJ 
and in Figure [S] we plot the results of the bias fitting. 

Rigorously, Equations [3] and f3] (and [TTJ and I12p only 
hold for a linear matter po wer spectrum (calculated using 
CAMB (jLewis et al.l 120000 ). as we have assumed the red- 
shift dependent power spectrum was separable: P(k, z) = 
D 2 (z)Pu„(k). We observe that using a linear power spec- 
trum, one can fit the ACF well for linear scales, i.e. for 
£ < 30. How ever we find t hat u sing a non-linear power spec- 
trum (from ISmith et ail (|2003l ). as implemented in CAMB) 
provides a better fit to the data, and to higher I, suggesting 
that using P(k,z) = D 2 (z)P non -ii n (k) can be considered a 
valid approximation to about £ ~ 50; since the bias does 
not change on the depth of our survey, we do not expect its 



/sky 21 + 1 



c, 



(14) 



and d = (C gg — C gg ); / s k y is the fraction of sky observed; C gg 
is our observed overdensity ACF after shot noise substrac- 
tion, given by 1/N where N is the mean number of galaxies 
per steradian; C gg is the theory (Equation HI and I12|l . Note 
that the covariance matrix M, and its determinant, depend 
on 6 g . 

We fit for the bias using our fiducial model (flat universe 
with tt m = 0.30, Q b = 0.05, h = 0.7, a 8 = 0.75) and we plot 
the resulting likelihood curve in Figure [6] The best fit value 
for the multipole range £ = 1-50 is b g = 1.40±0.03 (to 2 d.p) 
at la (fitting for a gaussian) . If we remove non-linear scales 
to the bias fitting, and consider only the multipole range 
t = 1 — 25, the error on the bias increases slightly while its 
value remains roughly the same (6 g ,im = 1-38 ± 0.05). We 
find our results in Section [7] do not differ significantly for 
these two values, and herein we take b g — 1.40 ± 0.03. 

The determination of the bias is particularly sensitive 
to as because they both act as overall normalization factors, 
and in linear theory C gg oc (b g as) 2 - 
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Figure 6. The (unnormalized) likelihood of b g from fitting to 
the 2MASS galaxy auto-correlation function for our flat fiducial 
cosmology (f2 m = 0.30, Q b = 0.05, h = 0.7, a 8 = 0.75). The best 
fit bias is f> g = 1.40 ± 0.03 from fitting for the multipole range 
£ = 1 — 50. When fitting over a strict linear range (I = 1 — 25) the 
best fit bias changes slightly to b g ,ii n = 1-38 ± 0.05 which does 
not affect the x 2 values reported in Table [2] The auto-correlation 
C gg is shown in Figure[5]with the fit b s = 1.40 . 



6 CROSS-CORRELATION METHOD 

We perform the cross-correlation in harmonic space. We 
have used masks on all the maps, and thus we cannot obtain 
true values of the multipole coefficients ae m , as the power 
of the harmonic estimator will be reduced and correlations 
will be induced between multipoles. However, we scale for 
the loss of sky cover and use a full covariance matrix to ac- 
count for the correlation between bins. We have obtained 
the spherical harmonic coefficients of our four 2MASS K20 
shells, and our four WMAP3 maps, as outlined in Section [5] 
We perform the cross-correlation: 

C MQ = irrr~~r\ V Re(a 3 lm aJ; n ) (15) 



(2£ + 1) ^ 



We further bin the data, using logarithmically spaced 
bins and I ^ 3. We avoid £ — 2 due to its anomalously low 
power in the CMB. For the analysis we use 5 bins 3 £ ^ 30, 
as this is where the ISW signal is expected to dominate. In 
Figures [7] and [8] we plot the correlation using 6 bins with 
3 < I < 200. 

In Figure[7]we compare the C 9 t results of the WMAP1 
and WMAP3 V-band data. Surprisingly we see a slight 
change, especially in the first point which corresponds to 
£ — 3 — 5. The power in the WMAP 3 maps has changed ve ry 
little (see for example Figure 19 in Irlinshaw et a l. (2006)), 
and thus the difference must be due to a slight change in the 
structure of these multipoles, p erhaps due to the improved 
gain model l|Jarosik et al.ll2006l ). Note that although the re- 
ported power spectrum has changed at \ow-£ this is due to 
change in the likeliho od analysis, rather than a change in 
the underlying data (|Hinshaw et aljfeood ). 

We will be comparing two hypotheses: a null hypothe- 
sis of no cross-correlation, and that returned by Equation [3] 
(and llip for our fiducial ACDM cosmology. To assess the fit 
we will compare the x 2 values returned by the two hypothe- 
ses. Similarly we can compare the unnormalized evidence £, 
where — 21n£ = \ 2 (if au parameters are fixed), or given by 
Equation [18] (if we marginalize over a parameter, where b s 
in the equation can be replaced by any parameter) . For this 
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Figure 7. The cross-correlation C 9 t (a*K) of galaxies from the 
2MASS catalogue, and the CMB from the V-band WMAP1 (open 
squares) and WMAP3 (filled squares) results. See Figure [8] for 
error bars. The power in the WMAP3 maps have changed very 
little (see for example Figure 19 in Hinshaw 2006), yet we see a 
slight change in the cross-correlations, especially in the first point 
which corresponds to £ = 3 — 5. The difference may be due to a 
slight change in the structure of these multipoles, perhaps due to 
the improved gain model (Jarosik 2006). 



we use a covariance matrix, estimated from simulations, and 
thus we account for the correlation between £ bins and those 
between K- shells. 

We make 500 simulations of a Gaussian CMB using 
the best fit theoretical Ctt from WMAP3. We apply the 
Kp2 mask, and correlate them with the four 7f-shell 2MASS 
maps. By not varying the 2MASS maps we have slightly un- 
derestimated the errors, because we have not accounted for 
the cosmic variance of the 2MASS data. We also note that 
the ISW signal in not inbuilt in the simulations; this should 
be inconsequential since this signal is negligible compared 
to the cosmic variance of the CMB. 

We compute the x 2 to find the favoured model and 
the improvement A(— 21n£) = A(x 2 )- As usual we have 
X 2 = d T Af _1 d, d = (C g T — C 9 t), where C g r is calculated 
for the two models, one with no correlation, i.e. C g r = and 
another with correlation due to ISW effect, given by equa- 
tion [3] (and llll) for our fiducial cosmology. The covariance 
matrix is defined as 



My = ((d i -(d i »(d j -(d j ))) 



(16) 



and is calculated from simulations in which the galaxies and 
the CMB are uncorrelated. We use 5 logarithmically spaced 
bins at low-i! where the expected signal dominates (i = 3 — 
30), and we include all four K20 shells in the analysis, thus 
i,j — 1,..,20. The fact that cosmic variance of 2MASS is 
not included suggests our simulations will underestimate the 
errors. We perform a consistency check on our covariance 
matrix by analytically estimating the diagonal elements of 
the covariance matrix: 
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Figure 8. Cross-correlation C 9 t(pK) results for the ILC (small 
triangle) , Q (open triangle), V (open square), and W (open pen- 
tagon) CMB maps with different magnitude bins of the 2MASS 
galaxy surveys. The dashed lines are la error bars about the null 
hypothesis, as evaluated from simulations. An ISW effect is ex- 
pected to be achromatic, which is what we observe, but the null 
hypothesis is not ruled out. 



\C gT ) 



f aky (2e + i) 



+ C„„Ctt) 



(17) 



which is the general form of Equation [14] We find that error 
bars calculated using Equation [17] are larger but of the same 
order of magnitude than those estimated from simulations. 



7 CROSS-CORRELATION RESULTS 

In this section we discuss the significance of the cross- 
correlation results and determine an upper limit on £Ia- 



7.1 Null Hypothesis 

In Figure [8] we plot the results of the cross-correlation, for 
6 logarithmically separated bins between f=3-200, with la 
error bars. The cross-correlation is achromatic, indicative of 
an ISW type cross-correlation. However the results appear 
completely compatible with the null hypothesis, and in fact 
scatter much less than expected. 

We use five logarithmically separated bins between £—3- 
30 to compute the — 2 In £ values for the (ILC,Q,V,W) maps. 
These values can be found in Table [2] (model 1) . By com- 
paring to simulations, we find these values are low at the 
~ 95% level. 



7.2 ACDM Fiducial Model 

To rule out the null hypothesis, at face value without consid- 
ering a competing theory, we would actually need a high \ 2 
value. However, we are comparing two theories and thus can 
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Figure 9. The observed C g T (a*K) (squares) from the V-band 
CMB map and the furthest magnitude shells of the 2MASS Sur- 
vey (13 < K20 < 13.5 and 13.5 < K20 < 14), with la cosmic 
variance. These two shells contain the most galaxies. The solid 
line represents the theory for our fiducial model (flat universe 
with Q m = 0.30, Q b = 0.05, h = 0.7, <r 8 = 0.75) and b g = 1.40. 
We have used 5 bins between I = 3 — 30 which corresponds to 
scales for which the ISW signal dominates. 



ask the more subtle question of which the data prefer. To do 
this we can consider the ratio of the evidences, A( — 21n£'). 

In Figure[9]we plot the results again for the two furthest 
redshift shells, which contain the most galaxies, for 5 bins 
at low-^ where the expected signal dominates. We also show 
the theoretical signal expected from our fiducial ACDM. 

Using these bins and the four K20 shells, the — 21n£ 
values found for the ISW theory, using b g — 1.40 and our 
fixed fiducial cosmology are reported in Table [2] (model 2). 
There is an improvement of A(— 21n£) ~ 1.5, and thus we 
confirm that the data prefer the ISW theory to the null hy- 
pothesis (of no correlation, or equivalently Einstein-de Sit- 
ter). However rule of thumb has it that 'strong' evidence is 
A(— In £ ) ^ 3, and thus this improvement is not compelling. 

We widen the uncertainty around the bias value, using 
the error bars from Section[S] There we found P(b s ) was ap- 
proximately a gaussian with ((J,, a) = (1.40,0.03), although 
we use the curve from Figure [6] for our marginalization: 



£= I P(data|theory,&g)P(6 g )db g 



(18) 



where we have assumed uniform priors on the theory. In 
Table [5] we record our — 21n£" results for a ACDM with 
a prior on the bias (models 3), which are identical at 2 
d.p. with those obtained from fixing b g — 1.40. We get 
very similar results if we use the 2D probability distribu- 
tion function P(b s , as) and marginalize over both b g and erg 
(using a gaussian with \x a i = 0.56 and a a i — 0.08). The 
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Table 2. Log evidence (— 21n£) values for cross-correlation of each WMAP3 maps (V, W, Q and ILC) with four 2MASS magnitude 
shells, using different model assumptions and priors. 
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Null Hypothesis 


ACDM 


ACDM 


ACDM 


ACDM 




of no correlation 


b g = 1.40 


prior on b g 


Q A = 0.85 


marginalized 












over Qa 


ILC 


11.3 


9.7 


9.7 


7.3 


10.1 


Q 


12.1 


10.4 


10.4 


8.1 


10.9 


V 


11.0 


9.5 


9.5 


7.4 


10.0 


w 


10.8 


9.1 


9.1 


6.9 


9.6 



Log evidence (— 21n£) values for cross-correlation of each WMAP3 maps (V, W, Q and ILC) with four 2MASS magnitude shells, using 
different model assumptions and priors. In all cases trg is taken to be 0.75. Model 1 is the null hypothesis of no correlation. In model 
2, we have considered a flat universe with Q m = 0.30, C;, = 0.05, h = 0.7 and b g = 1.40 (for b g = 1.38 only the W map result changes 
from 9.1 to 9.2). In model 3, we have widened the prior used on fe g to that which we obtained by fitting the 2MASS ACF on scales 
1 = 1 — 50, i.e. b g = 1.40 ± 0.03 (the results are unchanged when using b g .n n = 1.38 ± 0.05, the bias obtained from linear scales only, 
i.e., 1 = 1 — 25). Model 4 is best fit Qa = 0.85 in a universe with no curvature. In model 5, we have marginalized over Ha using the 
likelihood in Figure 1101 assuming a flat geometry and a uniform prior on Qa in the range [0;0.95]. The evidence ratio for model 1 (null 
hypothesis) and 5 (marginlized ACDM) is A(— 21n£) rs 1.1 which means the data prefer a ACDM cosmology to the null hypothesis, 
but only marginally. For models 1, 2 and 4 the evidence is related to the x 2 by — 21n£ = x 2 ■ 



results are also unchanged if we use the strictly linear bias 
6 g ,lm = 1.38 ±0.05. 



7.3 Assessing the goodness of fit 

By comparing our \ 2 values to those obtained from simu- 
lations, we find that they are low to ~ 95%, which is also 
evident from the lack of scatter in Figure [8] However, we are 
interested in maximising the evidence, or equivalently min- 
imising the x 2 1 an d thus we find that moving from the null 
hypothesis to the ACDM model, the fit is improved (raising 
the interesting question: is minimising the x 2 always appro- 
priate?). 

In the above analysis, we consider 20 correlated data 
points (five angular points in each of the four radial shells), 
and calculate the exact x 2 using a full covariance matrix. 
We can also consider each magnitude shell separately and 
calculate the \ 2 f° r each one, as a consistency check. The 
X 2 values obtained for the null hypothesis are (1.70, 0.67, 

0. 94, 2.36) going from the nearest to the furthest K20 shell, 
and for our fiducial model we find (1.72, 0.58, 0.83, 2.16). 
When the data are thus considered, our fiducial model is not 
always a better fit than the null hypothesis, and data in the 
closest shell prefer the null hypothesis. 

7.4 Upper Limit on Q, A and Marginalization 

Above we have compared our fiducial cosmology to the null 
hypothesis. Alternatively we can chose to consider only our 
fiducial ACDM model and use the cross-correlation to con- 
strain its parameters. We vary Qa, and Q m keeping all other 
parameters fixed (f2 m = 1 — Qa), and fit it to the measured 
correlation. 

In Figure [10] we plot the resulting likelihood. We find a 
best fit of Qa = 0.85, and upper limits of 0.87, 0.89, 0.90 at 

1, 2, 3(J respectively. Corresponding — 21n£ values can be 
found in Table [2] (model 4) . 

This relatively high value for Qa is in good agreement 
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Figure 10. The (unnormalized) likelihood of Ha from fitting to 
the observed C g T results (V band). The maximum is at Qa = 
0.85 Horizontal lines indicate the 1,2,3 a levels, equivalent to the 
decrease in the likelihood value. We find the corresponding upper 
limits on Ha at 0.87, 0.89,0.90. Corresponding x 2 values can be 
found in Table [2] 



with other studies of the ISW effect: ICabre et all (|2006f ) find 
Qa — 0.8 — 0.85 at the la level. As can be seen, the null 
hypothesis (Qa = 0) is less than 2a away from the best fit 
result, confirming that we cannot confidently rule it out. 

In the last column of Table [2j we present the marginal- 
ized evidence for our ACDM model (model 5), where the 
prior on the theory is flat over Qa — — 0.95, and the like- 
lihood of our model is taken from Figure [lOl For this model 
A(— 21n£) ~ 1.1 on average, so the data prefers a ACDM 
model to the null hypothesis, but only marginally. 

Bearing in mind this issue, and the fact that the % 2 
values are low for all models, any claim for an ISW detection 
using 2MASS and WMAP3 remains tentative, and indeed 
our results are also consistent with the null hypothesis of no 
correlation within the 2a level. 
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8 AXIS OF EVIL 

We briefly consider other possible sources and statistics of a 
positive cross-correlation between the CMB and LSS. As dis- 
cussed in Section [1] interesting anomalies dubbed the "Axis 
of Evil" (AoE) have been observed in the WMAP data of 
the CMB that indicate a possible departure from statistical 
isotropy (|Land fc Mague iio 2005a). Here we give a brief de- 
scription of the AoE anom aly, but the reader should refer to 
lLand fc Magueiio! (|2005al ) for a detailed description. 

In harmonic space, a statistically isotropic map is ex- 
pected to have independent I modes. For a given multipole 
I, the power should also be distributed randomly amongst 
each m-mode, but it is possible to rotate any map so that 
the power in a given multipole is mostly in a given m-mode. 
However, if for a given frame, several multipoles have power 
mostly distributed in one m-mode, the n the map can no 
longer be considered statiscally isotropic. lLand fc Magueiiol 
|2005al ) found that there existed a set of nearly identical 
frames (with axis dubbed the "Axis of Evil") in which the 
multipoles of the WMAP1 maps showed phase correlations. 

Possible sources of these anomalies are: foreground con- 
tamination; astrophysical effects (e.g., lensing and mov- 
ing cluster effect); alternative cosmological paradigms. If 
the AoE was in some way due to local inhomogeneities 
(e. a.. rValel (|2005l )l then one might also expect a C g r cross- 
correlation on these large scales. We investigate the AoE 
anomaly in two ways. Firstly, we rotate the 2MASS data to 
the frame of the WMAP AoE frame and search for anoma- 
lous phase correlations in 2MASS in that frame. Secondly, 
we rotate the 2MASS data in all possible directions, search- 
ing for an AoE type anomaly in all frames. 

For the first test, we examine if the AoE signal observed 
in the CMB is also present in the 2MASS data. In Figure [TT] 
we plot the power observed in each m-mode, for the AoE 
frames returned by the CMB. For each multipole, this is 
the frame where one m-mode dominates. We plot the power 
ratio 2\ai m \ 2 /Ce (without the 2 factor for m = 0) evaluated 
in the AoE frames, for the 2MASS (12 < K 2 p < 14) as well 
as the CMB (using the cleaned ILC map of iTegmark et ail 
(2003)). By definition, in these frames the CMB observes 
a pattern of one m significantly dominating each multipole. 
We do not observe a similar pattern in the 2MASS data, and 
thus conclude that if there exists a source responsible for the 
AoE features, then is is unlikely to be the same source of the 
C 3 t correlation observed above. 

For the second test, we looked at 2MASS indepen- 
dently over I — 2 — 20 and found no evidence for any AoE 
type structure (a correlation of multipole frames as defined 
above). Clearly the 2MASS catalogue is highly anisotropic, 
with structures such as the Supergalactic Plane visible by 
eye as well as the Galactic mask. Is it not strange that a 
statistic that measures evidence for deviations from statis- 
tical isotropy in the CMB does not find any such feature 
in 2MASS? The AoE statistic highlights a very particu- 
lar type of phase correlation, in harmonic space, and is by 
no means a conclusive test for general deviations from sta- 
tistical isotropy. The fact that a clearly non-gaussian and 
anisotropic map, such as 2MASS, does not return a anoma- 
lous AoE signal for I = 2 — 20 perhaps highlights the weak- 
ness of this statistic, and throws caution at how one de- 
fines and selects appropriate statistics. Perhaps the non- 
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Figure 11. The fraction of power in each m-mode for multipoles 
£ = 2,3, 4, 5. The ae m s are computed in the Axis of Evil (AoE) 
frames returned by the CMB. By definition the multipoles of the 
CMB are dominated by one m-mode in these frames (shaded). 
We do not observe a similar pattern in the 2MASS results (solid 
line). Negative values of m are not shown as a^_ m = aj m . 

gaussianity is 'washed out' by the depth of the survey - 
anomalies might be more significant in the shallowest part 
of the survey. However, there is currently no conclusive way 
to test a map for deviations from statistical isotropy. 



9 DISCUSSION 

We calculate the cross-correlation between the 2MASS 
ga laxy survey and the W MAP3 data. This updates the work 
of lAfshordi et alj (|2004h who cross-correlated 2MASS with 
WMAP1. 

The cross-correlation signal expected in a ACDM Uni- 
verse scales with the linear galaxy bias, b g , and in linear the- 
ory with the product (6 g crg). We fix as — 0.75, and use a flat 
fiducial cosmology (based on inflation). We use fl m = 0.30 
(based on supernovae la data), Qt = 0.05 (based on Big Bang 
nucleosynthesis), h — 0.7 (based on the HST key project), 
n = 1 (i.e., the Harrison-Zeldovich spectral index). 

Fitting the fiducial cosmology to the angular auto- 
correlation function of the 2MASS galaxy survey for 12 < 
K20 < 14, yields a linear bias value of 6 g = 1.40 ± 0.03, as- 
sessed from multipole scales i = 1 — 50 (6 g ,ii n = 1-38 ± 0.05 
assessed from multipole scales I = 1 — 25). As the 2MASS 
galaxy survey is shallow, we assume b s is constant with red- 
shift over the depth of our survey. 

The measured cross-correlations obtained from four dif- 
ferent WMAP maps (V, W, Q, and ILC) and four dif- 
ferent magnitude shells of 2MASS (12 < K 2 o < 12.5, 
12.5 < K20 < 13, 13 < K20 < 13.5, and 13.5 < K 20 < 14) 
show an achromatic signal, as expected from an ISW effect. 
However, the observed signal is also within the la error bars 
obtained from cross-correlating random simulations of the 
CMB with 2MASS data. This means the data are consis- 
tent with the null hypothesis of no correlation. 

We compare our observation with the ISW signal ex- 
pected in our fiducial ACDM model, with a fixed bias and 
an uncertaintly around the bias value which we found when 
fitting the auto-correlation function. The ACDM model finds 
a lower chi-squared, and thus is a better fit compared to the 
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null hypothesis, but in both cases of treating the bias the 
evidence change is only A(— 21n£) w 1.5. Whichever model 
is considered - the null hypothesis or a ACDM universe - the 
X 2 values are low when compared to simulations. 

Varying the Dark Energy density component of the 
ACDM model and assuming flatness, we find that the data 
prefer a high value of Q A = 0.85, with fi A < 0.89 (95% CL). 

Thus, there is a higher correlation between WMAP3 
and 2MASS than that expected by a ACDM universe with 
SIa = 0.7. However, the observed cross-correlation may not 
only be due to an 1SW effect. Other signals might contribute 
to the correlation: positive curvature (even if very small) or 
cosmic magnification (higher redshift galaxies can contribute 
more than expected if they are lensed by the lower redshift 
galaxies and their luminosity function is such that the effect 
leads to a positive correlation). 

When we marginalize over Qa, the evidence ratio be- 
tween the null hypothesis and ACDM becomes A(— 2 ln£) ~ 
1.1 so we can say that the data prefer a ACDM universe, but 
only marginally. In any case, the correlation observed is con- 
sistent with both hypotheses within 2a. 

We also investigate "Axis of Evil" (AoE) type anoma- 
lies, which detect phase correlations between different mul- 
tipoles in harmonic space. These phase c orrelations are not 
expect ed in a statistically isotropic map (|Land fc Magueiiol 
|2005a|). We do not observe an AoE type of structure in the 
2MASS catalogue. As non-Gaussian features are expected 
in the LSS we feel this result raises issues about the use of 
the AoE statistic as a general test for statistic anisotropy. 

We do not observe correlations between the CMB frac- 
tional power distribution as measured in its AoE frame and 
that of 2MASS, constraining the possible explanations of 
the low £ anomalies in the CMB. 

Future spectroscopic and photometric galaxy redshift 
surveys (e.g., the Dark Energy Survey, the Wide Field Multi- 
Object S pectrograph) wi ll yield more galaxies out to higher 
redshifts. lAfshordil (|2004 ) showed that an all sky survey with 
10 million galaxies and uniform sky coverage between < 
z < 1 would lead to a detection of the ISW effect at the 5a 
level. It remains to be assessed what combination of depth 
and sky coverage is optimal for detecting the ISW effect. 

We highlight that a claim for an ISW detection could 
be greatly weakened if one considers some of the uncer- 
tainty around the cosmological parameters. For our 2MASS- 
WMAP3 correlation the signal-to-noise ratio is poor, how- 
ever with a stronger data set a full Markov Chain Monte 
Carlo exploration of parameter space should be done. 
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